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A spatially one dimensional coupled map lattice possessing the same symmetries as the Miller 
Huse model is introduced. Our model is studied analytically by means of a formal perturbation 
expansion which uses weak coupling and the vicinity to a symmetry breaking bifurcation point. In 
parameter space four phases with different ergodic behaviour are observed. Although the coupling 
in the map lattice is diffusive, antiferromagnetic ordering is predominant. Via coarse graining the 
deterministic model is mapped to a master equation which establishes an equivalence between our 
system and a kinetic Ising model. Such an approach sheds some light on the dependence of the 
transient behaviour on the system size and the nature of the phase transitions. 



I. INTRODUCTION 



Since the middle of the seventies the investigation of deterministic chaos has become one of the prominent fields 
in science, especially physics. A lot of knowledge has been gained since that time, in particular for low degree of 
freedom systems [Q, and a whole machinery of tools has been developed for the diagnostics of chaotic motion. We just 
mention Lyapunov exponents and fractal dimensions as the most popular quantities. Parallel to these developments 
the question has been raised to which extent the number of degrees of freedom enters the business. Unfortunately, 
much less progress has been achieved in this direction. Only few results are available and most of them are bound 
to the investigation of model systems. Within that context coupled map lattices (CMLs) have been introduced at 
the end of the eighties as a widely studied model class fn such models local chaos is generated by a chaotic 

map which is placed at each site of a simple lattice. Spatial aspects are introduced by coupling these local units and 
special emphasis is on the limit of large lattice size where the dynamics becomes high dimensional. 

There is just one class of many degree of freedom systems which is fairly well understood, namely statistical 
mechanics at and near thermal equilibrium. Unfortunately, the systems studied in the field of space time chaos are 
often far from equilibrium so that the tools of equilibrium statistical mechanics may fail. Nevertheless, the reduction to 
relevant degrees of freedom, sometimes called coarse graining, may be equally successful in both areas. By elimination 
of irrelevant degrees of freedom one maps the microscopic deterministic equation of motion to a stochastic model where 
the noise captures the irrelevant information. Such a concept, well developed in equilibrium statistical mechanics, has 
also been used in nonlinear dynamical systems; introductions can be found on the textbook level 0]. In a rigorous 
approach coarse graining is performed by suitable partitions of the phase space and there are results for particular 
coupled map lattices available (cf. (^||). Unfortunately, such schemes are limited to some perturbative regime and 
are technically extremely difficult to apply. Henceforth, sometimes more physically motivated coarse grainings are 
used relaxing the amount of rigour a little bit. 

The just mentioned statistical methods become especially relevant in the study of phase transitions in CMLs [p|. 
Qualitative changes in the dynamical behaviour may be related to phase transition like scenarios in the corresponding 
coarse grained description. Prominent examples for such phenomena occur in the models introduced by Sakaguchi M 
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and Miller and Huse M] ■ To keep the paper self contained and as a motivation for the construction of our model we 
shortly review basic features of the latter model. 

In order to mimic a phase transition in a two dimensional Ising model, the chaotic antisymmetric map depicted in 
figure was placed onto a square lattice and coupled to its four nearest neighbours 



= (l-e)^y + - 



k,l=±l 



(1) 



Performing a coarse graining according to the sign of the phase space variables 
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(2) 



numerical simulations indicate a phase transition if the coupling strength exceeds a critical value e C rit. ~ 0.82 (cf. figure 
|l|). Extensive numerical simulations indicate |p"l|| that the phase transition is continuous. However, it is doubtful 
whether the transition belongs to the Ising universality class, because the results for the critical exponents are 
inconclusive. In particular, their values depend on whether the CML is updated synchronously or asynchronously. 
One can summarise that the phase transition of the Miller Huse model is still far from being understood, in particular 
since no quantitative description of the spin dynamics could be derived. In order to reach some progress in this 
direction we here introduce and investigate a slightly different model system with analytical methods. 




FIG. 1. The single site map <f) of the Miller Huse model. The inset shows a typical pattern for e > e cr it. (white/black 
a tJ = ±1). 

Section [l| introduces our model as well as the setup of the perturbation expansion. For the latter purpose transitions 
between sets of a suitable partition are defined. These transitions are studied in detail in section III. Herewith, the 



bifurcation diagram of our model will be developed in section IV and analytical expressions for the bifurcation lines 
are calculated in perturbation theory. Section [V] is devoted to a systematic coarse graining of the dynamics on the 
basis of the just mentioned partition. On that level the dynamics is described in terms of a master equation which 
corresponds to a particular class of kinetic Ising models. It constitutes the basis for the investigation of the transient 
behaviour in section VI. Finally, the main results of this work are summarised. The appendices are concerned with 
parts of the perturbation expansion, but more details can be found in |12| . 



II. THE MODEL 

Let us first consider the single site map. It consists of a deformed antisymmetric tent map fg, which is linear on 
three subintervals of [—1, 1] 

!— 2 — x/a , if x E[— 1,— a] ^ 
x/a , if x€(-a,a) , a := - — - . (3) 
2 - x/a , if x S [a, 1] 10 

Because of = S the parameter 6 determines whether transitions between the intervals [—1,0] =: J(— 1) and 

[0, 1] =: J(+l) are possible. Note that the Miller Huse map is obtained as a special case, cj) — fs=-i- The introduction 
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of a in eq. (|3|) ensures that the modulus of the derivative of fs is constant on the whole interval. Figure g shows the 
function fs for small positive and negative S. For 5 > the single site map has two coexisting attractors, the intervals 
[— 1. S] and [5, 1], whereas for S < only one attractor, the interval [—1, 1], is present. 




FIG. 2. The deformed antisymmetric tent map fs. 



The CML which is studied in this article is defined on a one dimensional lattice (chain) of length N. Nearest 
neighbours are coupled in a standard "diffusive" way with periodic boundary conditions 

T e , 5 : [-1, [-!,+!]", 

[T e , 5 (x)] 4 := (1 - e) fs( Xi ) + | (/,(a;,_i) + fs(x i+ i)) . (4) 

The parameter e denotes the coupling strength. Because of the single site map and the diffusive coupling the CML 
T e .,5 has the symmetry T e x) = — T C! a(x). Furthermore, translation invariance on the one dimensional lattice 
holds, because periodic boundary conditions have been imposed. 

Since we are going to perform a perturbation theory with e, \S\ <C 1, we first consider the CML with e = 8 = 0. 
In this case the model can be solved trivially. The non-deformed antisymmetric tent map /o has the two attractors 
,/(— 1) = [—1, 0] and = [0, 1]. Therefore, N uncoupled maps fo have 2 N coexisting attractors, each one an N 

dimensional cube of edge length one 

I a := J(ai) x J(a 2 ) x . . . x J(ccn) ■ (5) 

We distinguish these cubes I a by an N dimensional index vector a. = [otx, a>2, ■ ■ ■ ajv) where an € { — 1, +1}- The 
natural measure on each cube is the Lebesgue measure. As we will see these cubes become important building blocks 
of the perturbation theory and the starting point of a coarse grained description of the CML T e j. 

From a dynamical systems point of view we are mainly interested in ergodic properties of the CML, i. e. the number 
of coexisting attractors and their location for given small parameters e, 5. An important observation is that in the 
perturbative regime a typical orbit stays for many iterations within a cube I a , before it possibly enters another cube 
Ip. Therefore, in perturbation theory any attractor of the CML T Cj( 5 is a union of cubes I a , if one neglects sets with 
volume 0(e, S). Hence, the dynamics is sufficiently characterised by transitions I a — » Ip between cubes. 

Of course we have to be more definite with what we mean by a transition. In order that a phase space point can be 
mapped from a cube I a to a cube Ip (a ^ (3) the image of the former has to intersect the latter. Hence the overlap 
set 

OV a ,p:=T etS (I a )nIp (6) 

plays an important role. A necessary condition for a point to migrate from I a to Ip is a non-empty overlap set 
OV at p. Since in perturbation theory the set T t ^{I a ) is a weakly deformed cube I a , the set OV a ,p can at most have 
a volume of size C(e, S). However, the condition on the overlap set is far from being sufficient because one has to 
ensure that typical orbits can reach this set upon their itinerary. For that purpose two additional conditions have to 
be imposed. 

First, we have to ensure that points from the inner partj^Jof I a reach the overlap set. For that reason we consider 
the pre-images of OV a . p of various generation that are contained in I a 



For our perturbative treatment we define the inner part as the set of all x G I a which have at least a small fixed positive 
distance d from the boundary, where the quantity d does not depend on the expansion parameters e and 5. 
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T-}(OV a ^) := {xe/ a T £ ,*(x) G OV a ,p} 

T- 5 fc (01/ a ^):={xe/ a |T £ , 5 (x)GTj fc - 1) (OK l , / 3)}, k = 2, 3, . . . . (7) 

For some finite k the pre-image set (OV^^g) should intersect the inner part of I a , so that points from the inner 
part of I a can reach the overlap set OV a .p^\ 

The points of the set OV a: p are near the surface of the cube Ip within a distance of order 0(e, S). The second 
condition demands that points from a subset of OV aj p with finite Lebesgue measure reach the inner part of the cube 
Ip directly under further iteration. The two conditions for a transition — * Ip ensure that the transition is possible 
for a set of finite Lebesgue measure that is located in the inner part of I a . 



III. TRANSITIONS IN PERTURBATION THEORY 



In what follows we consider the CML T e _$ for arbitrary but fixed lattice size N. We would like to know which 
transitions — > Ip are possible for given parameters e, S. In the spirit of perturbation theory we confine ourselves 
to dominant transitions. Those are transitions where the cubes I a and Ip share an (N — 1) dimensional surface. 
Then, the volume of the overlap set OV at p can be greater by a factor 1/e or 1/|<5| in comparison to the case without a 
common surface. Consequently, the N dimensional index vectors a and (3 only differ in one component, the transition 
index on. In such a transition I a — > Ip the Xi coordinate of the phase space orbit {x 4 } changes its sign. Transitions of 
higher order in which two or more coordinates simultaneously change their sign will not be considered in this article, 
because their rates are smaller by a factor of the order 0(e, S) in comparison to the dominant transitions. 

In perturbation theory, for a dominant transition only the neighbouring indices of the transition index, ai-\ and 
Of +i, are relevant, because of the nearest neighbour interaction of the map T e ,6 (cf. eq. (^)). In addition, the 
influence of the two neighbouring coordinates and Xi + \ on the Xi coordinate is predominant for a finite number 
of iterations, since interactions with lattice sites farther away are suppressed by the small coupling strength e. More 
precisely, within first order perturbation theory the overlap sets OV a ^ p and their pre-image sets can be approximated 
by the following product sets (cf. appendix |A|) 

(IV a - OV {3) x r^- 3 ) 

<~>Va,p — ^ v a i _ 1 a i a i j rl ,l3i^ 1 flil3i +1 A J ai o 2 ...a,_ 2 a i+2 ...aiv ' 

T; k s (OV a ,p)= [T[ 3 i]^(0^ iQiai+ii/3i _ i/3i/3i+i ) x l£-*U_ 2ai+2 ... aN , k>l . (8) 

(3) 

Here OV^ i _ l aj 1 /3 ._ 1 p. p. 1 denotes a three dimensional projection of the full overlap set which contains the coordi- 

nates Xi_i, Xi and Xi+i, and T^ { denotes the map lattice for N — 3. The (TV — 3) remaining coordinates are contained 

in the (N — 3) dimensional cube I^a2---ai- 2 a i+2 ...a N ■ Effectively, we have herewith reduced the transition in a map 
lattice of size N to a transition in a map lattice of size three, because (N — 3) coordinates play only a spectator role. 
Put differently, the CML T £i a reaches already its full complexity for N = 3, if one stays in the perturbative regime. 
For symmetry reasons one can identify three different types of transitions I a — > Ip: 

Type (a): the three indices oti—i, on and ai+i are equal, e. g. 

+1,4-1,4-1,... — * -1,+1,... ■ 

Type (b): the two neighbouring indices and oii+i are different from each other, e. g. 

I..., — i,+i,+i,... — » I...,-x,-x,+\,... 

Type (c): the neighbouring indices aj_i and a^+i differ from the transition index o^, e. g. 

' .11-1 — > ' .1.1.1 



2 Since for e = 5 — the natural measure on each cube is the Lebesgue measure, in the perturbative regime the map T e ,s 
distributes the points of an orbit rather uniformly within a cube I a . Therefore, in determining the orbit dynamics it suffices 
to use topological methods like the calculation of pre-image sets. 
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Transitions of type (c) are inverse to those of type (a) . 

Because of the conditions mentioned in the last section transitions are possible only if the deformation is small 
enough, 8 < <S cr ,*t. (e). Within perturbation theory we obtain for the different critical values 

2e 4e 
type (a) : S a = , type (b) : 6 b = — — , type (c) : S c = -— . (9) 

One might wonder why transitions (b) and (c) do not appear for negative 5 above the critical value. The main reason 
is that despite of the existence of a non-empty overlap set trajectories do not reach this overlap since there exists a 
forbidden region in phase space called the " blind volume" . Points belonging to the blind volume have no pre-images 
themselves. The blind volume is non-empty, since the map T e ,s is not surjective for finite coupling e. The actual 
calculation of critical 6 values necessitates rather involved geometric constructions in phase space, since one must 
determine the location of the pre-image sets TJg(OVa, p) in la- Hence, details are postponed to appendix [§]. The 
smaller the deformation parameter <5, the more transitions become possible as can be guessed from the geometry of 
the map (cf. figure On the other hand, increasing coupling constant e inhibits transitions, since eventually only 
the transition of type (a) remains feasible for fixed negative 8. Such an observation contradicts somehow the intuitive 
reasoning about a "coupling" of lattice sites. The inhibition effect for transitions is caused by the existence of a "blind 
volume" in the cube I a that grows with e (cf . appendix |b|) . 

At this stage some remarks on the accuracy of our perturbative approach seem to be in order. Since we neglect 
transitions of higher order our arguments are not rigorous. In fact for a real proof the complete absence of such 
transitions must be shown. For the case of two coupled maps, N = 2, such a step can be easily supplemented 
(cf. appendix |c|) and we infer that one might be able to perform similar but more involved computations in higher 
dimensional cases too. Nevertheless, even if these transitions are mathematically possible their effect may be small 
e.g. taking a time scale argument into account. 



IV. THE BIFURCATION SCENARIO 

Eq. @ determines four regions in the (e, 6) parameter plane where different transitions are possible (cf. figure^). 
Crossing these lines a bifurcation occurs. Between different regions the number of coexisting attractors and their 
location change. Determining attractors in the strict mathematical sense just from the knowledge of the dominant 
transitions faces however some problems. For, neglecting sets with volume C(e, 5) a union of cubes A, which an orbit 
can not leave through a dominant transition, is a candidate for an attractor. But, possibly an orbit can escape from 
this set through a transition of higher order in perturbation theory. Then the set A would not be an attractor in a 
strict mathematical sense. However, one can put forward the following time scale argument: in perturbation theory 
transitions through which an orbit can leave A occur on a rather large time scale in comparison to the relatively fast 
dominant transitions through which the orbit is pulled back to the set A again. Because of this intermittent dynamics 
the set A is a core region of a possibly bigger attractor, i. e. the sets A are the carriers of most of the natural measure 
of these attractors. For brevity we will call these sets A "attractors" in the following. 

If we neglect sets with volume 0(e, S) we can identify an attractor A with a union of cubes I a . 
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FIG. 3. Diagrammatic view of the bifurcation diagram for the CML T e ,s according to eq. (J9|). The numbers in the four 
parameter regions refer to the text. Gray-shading indicates the type of coupling in the corresponding kinetic Ising model, 
antiferromagnetic (light) resp. ferromagnetic (dark) (cf. section |v|). 

Region 1 (<5 > 0) : No transition I a — > Ip is possible. Therefore, each cube I a is an attractor so that there are 2 N 
coexisting attractors. 
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Region 2 (— 2e/3 < 5 < 0): Only transition type (a) is allowed. Hence, cubes I a are attractors such that a does 
not contain three successive or " — 1". With a combinatorial argumentation one can show for that for long 

chains (N ^> 1) the number of coexisting attractors increases like ((1 + V5)/2) N . 

Region 3 (— 4e/3 < 5 < — 2e/3): To determine the attractors in this region, it seems necessary to anticipate the 
coarse graining of the CML T C: ,5 which will be discussed systematically in section |v|. Analogous to eq. (|^) wc 
can view the index vector a of a cube I a as a spin chain of length N, where +1 and —1 are the possible spin 
states on each lattice site. In this way the three transition types (a), (b) and (c) translate into three different 
kinds of spin flips. For each spin chain one can define defects in the same way as in the antiferromagnetic Ising 
model. A defect ("1") occurs, if two neighbouring spins are aligned, and no defect is present, if the spins point 
in opposite directions. Then, the just mentioned spin flips translate into a dynamics of defects. 

Type (a): two adjacent defects annihilate each other, e. g. 

spin chain a : . . . + 1, +1, +1, —»... + 1, —1, +1, . . . 

defects in a : 1, 1, — > 0, 0, 

Type (b): one defect diffuses to a neighbouring lattice site, e. g. 

spin chain a. : . . . + 1, +1, —1, —>... + 1, —1, —1,... 

defects in a : 1, 0, — > 0, 1, 

Type (c): two adjacent defects are generated simultaneously, e. g. 

spin chain a : . . . + 1, —1, +1, —>... + 1, +1, +1, . . . 

defects in a : 0, 0, — > 1, 1, 



For the determination of the attractors in the present parameter region we consider an orbit {x*} of the CML 
which performs successive transitions I a — ► I p. Each transition changes the corresponding spin chain a and 
its defects. Since transitions of type (c) are forbidden, defects can diffuse and annihilate in pairs only, and the 
number of defects decreases monotonically. 

If the size of the system N is even the chain contains an even number of defects. An orbit {x*} migrates 
between cubes I a , until all defects have annihilated each other. Then, the orbit can not execute any fur- 
ther transition of type (a) or (b). Therefore, there are two attractors, the cubes 

and x,+i, Extensive numerical simulations indicate that for N even each cube 

I( +1,-1, +1,-1, +1,-1) an d ^(-i,+i,-i,+i,...,+i,-i,+i) constitutes an attractor in the strict sense, i. e. for 
S > — 4e/3 no additional transition of higher order perturbation theory is present (cf. also appendix 

For N odd the number of defects in a is odd. Consequently, at the end of the transient dynamics one defect 
remains. Since the defect can change its location via a transition of type (b), the attractor is the union of all 
2N cubes I a for which a contains a single " + 1 + 1" or " — 1 — 1" sequence. 

Since in both cases the ratio of the volume of the attractor to the volume of its basin of attraction becomes very 



small for A ^ 1 one expects long transients to occur. Section |V| is devoted to a more detailed study of the 
transient dynamics. Our argumentation has used the assumption that different transitions are not correlated. 
We will come back to this problem in the next section. 

Region 4 (S < — 4e/3): All three transition types are possible. Therefore an orbit {x*} can visit every cube I a , so 
that there emerges one attractor which encompasses all cubes. 



V. COARSE GRAINING OF THE CML 



Coarse graining the CML T ey s one passes from orbits {x*, t — 0, 1, 2 . . .} in phase space to symbol or spin chains 
{a*, t = 0, 1, 2 . . .}. The spin chain a* just indicates the cube which contains the phase space point x* at time t 
(cf. eq. (||) ) • If an orbit of the CML performs a transition I a — > Ip, the state of the spin chain changes from a. 
to f3. Since in perturbation theory an orbit typically circulates for many iterations within a cube I a , the sequence 
{a*, t = 0, 1,2...} has a constant value for long time before a spin flip occurs. Altogether, the CML is described by 
a stochastic spin dynamics. First we argue that the spin dynamics is Markovian for the following reasons: 
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• Two successive transitions I a — > Ip and Ip — > I~, are uncorrelated. In the perturbative regime an orbit performs 
a highly chaotic motion within the cube Ip for many iterations before the transition to the cube I-y occurs. 
Therefore, the memory of the preceding transition I a — > Ip is lost. 

• A transition I a — ► Ip is equally probable for each iteration step. For a transition the orbit point x* must 
hit a characteristic set in the inner part of the cube I a which consists of pre-images of the overlap OV a _p. 
During its stay within the cube I a the orbit is distributed uniformly within I a , since for e = 5 = the natural 
measure on I a is the Lebesgue measure. Consequently, the probability for the orbit to hit the characteristic set 
is independent of time. 

Since the spin dynamics is Markovian at least approximately, the probability p a (t) that the spin chain is in state a 
at time t obeys a master equation with transition probabilities w (/3 | a) for a spin flip a — > (3 

p a {t+l)=p a (t)+^[w(a\/3)-pf l (t)-w(p\a)'p a (j;)] . (10) 

In perturbation theory three types of spin flips occur. As already stated above only a single spin flips during the 
elementary process a — ► f3. From the study of the underlying CML T e j one can infer the following properties of the 
transition probabilities w (/3 | a). 

• Because of the nearest neighbour coupling in eq. (|]) and the direct product property (||) in the perturbative 
regime, the transition probabilities w(J3 \ a) depend on the three neighbouring spins only 

w(f3 | a) = w (3) (a,_i ^ a l+1 | a, a i+1 ) , (11) 

where on denotes the transition index. Hence the spin interaction is local. 

• According to eq. ( p"l| ) the transition probabilities w (/3 | a) do not depend on N. Therefore, they can be deter- 
mined for small systems, e. g. N = 3. 

• Transitions of the three different types have probabilities w a , Wb, and w c , respectively, which depend on the 
parameters e, 8. If 5 is greater than the corresponding critical value in eq. (^|), the respective transition prob- 
ability strictly vanishes. Lowering 5 the transition probability increases monotonously as can be shown by a 
rather subtle argument which uses the monotonous growth of the overlap set OV a .p. 

The dynamics resulting from the master equation ([Io|) is almost trivial in parameter regions 1 and 2, since at most 
the spin flip of type (a) is possible. Regions 3 and 4 are more interesting, because at least two different spin flips 
occur. Since we are mainly interested in large systems, we confine ourselves to N even in what follows. 

In region 3 spin flips of type (a) and (b) are possible. On the coarse grained level the attractors 
I( +1,-1,.. +1,-1 ) an d I(— i,+i,— i,+i,...,+i,— i,+i) are viewed as the two ground states of the antiferromag- 
netic Ising model. Hence, the ergodic dynamics of the CML T ei s corresponds to an antiferromagnetic Ising model at 
zero temperature. In parameter region 4 all three spin flips are possible. The (unique) stationary distribution of the 
master equation can be calculated with the ansatz that the weight of each state a solely depends on the number of 
defects. The result 

= ^exp L 8 J ^2aiOH +1 J (12) 

clearly can be cast into the form of a canonical distribution for a nearest neighbour coupled Ising chain. Here, c and 
Z denote the normalisation constants and for the temperature the relation 

follows. Taking J with modulus one the temperature depends on the ratio of the transition probabilities for generation 
and annihilation of two defects. It is finite throughout region 4. Ferromagnetic coupling, J = +1, is obtained for 
{w c /w a ) > 1, whereas in the opposite case {w c /w a ) < 1 antiferromagnetic coupling, J = —1, follows. Both cases are 
realized in the present parameter region, as can be seen in figure [| Here, the transition probabilities w a and w c were 
obtained numerically by analysing the transitions of a very long orbit of the CML with N — 3. Finally, it is easy to 
show that the stationary distribution ( |l2| ) of eq. ( |l0| ) obeys detailed balance 

w(/3 | a) ■ p% at = w{a | (3) ■ p s p tat , 

although the underlying CML describes a non-equilibrium process on the microscopic level. 
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VI. TRANSIENT DYNAMICS OF THE CML 



As mentioned above the transient dynamics of the CML is most interesting in the parameter region 3 for which the 

attractors It and 1/ +1,-1. +1) have large basins of attraction for large N. The 

mean transient time (T) can be determined by averaging the time until an orbit reaches one of the two attractors over 
many random initial conditions, i.e. initial conditions distributed according to the Lebesgue measure. The numerical 
simulation indicates a quadratic increase with the system size 

(T) (xN 2 (N > 1) . (14) 

Such a law can be understood from the coarse grained point of view. In parameter region 3 the spin flips of type 
(a) and (b) are allowed which cause the annihilation of two defects respectively the diffusion of one defect. Overall, 
the transient dynamics of the CML corresponds to a relaxation process towards one of the two ground states. Since 
diffusion is important for the relaxation the time scale grows with the second power of the length scale, i.e. the system 
size. 

To be more definite and in order to derive eq. ([l4]) formally we remind that the spin dynamics induced by the 
map lattice constitutes a kinetic Ising model with local spin flips. Models of these type have been introduced by 
R. J. Glauber in his celebrated article (l^Jl^]. So the coarse grained dynamics of the CML T t ,s belongs to a well- 
studied class of models. For the case of zero temperature, i. e. w c = 0, an exact analytical solution is available if an 
additional relation for the two remaining transition probabilities is imposed 

w a = 2wb . (15) 

Such a condition holds only on a subset of parameter region 3, namely on a line of e, S values. One can show with the 
help of results on the temporal evolution of correlation functions that the mean number of defects in the spin chain 
obeys 

1 N 

(# defects (i)) ~ -== -= , t > 1 . (16) 

y'STTWb Vt 

If one changes the system size from TV to k ■ N the time t has to be scaled by a factor k 2 in order to reach the same 
number of defects. Since the mean transient time (T) determines the scale for the annihilation of all defects, the scale 
argument implies relation (|l4|). 

When one relaxes condition ( |i"5| ) between transition probabilities the quadratic growth of the transient time with 
N in eq. ( [l4|) still holds as numerical simulations indicate. Such an observation is in accordance with the theory of 
dynamical critical phenomena |T5| ] . The latter implies universal scaling laws for relaxation phenomena at the critical 
point. At zero temperature the dynamics of a one dimensional kinetic Ising model is critical and the decay of defects 
is governed by the dynamical critical exponent z 

(#defects(*)>~^ . (17) 

z equals two for kinetic Ising models where the order parameter is not conserved |l5|,[l6| . Consequently eq. (|l4|) holds 
for a large set of parameter values in region 3. 



VII. SUMMARY 



We have introduced a coupled map lattice which was constructed in analogy to the Miller Huse model. By a 
perturbation expansion for weak coupling and in the vicinity of a symmetry breaking bifurcation of the single site 
map nontrivial dynamical behaviour has been investigated. Our approach was based on analysing geometric properties 
in phase space. Transitions between certain cubes which are the building blocks of a coarse grained description have 
been computed. A global bifurcation of the dynamics occurs if a transition becomes allowed or forbidden by a change 
of the parameters. Four parameter regions with different ergodic behaviour could be identified. As a surprising 
and counter intuitive feature of our map lattice we mention that increasing the spatial coupling inhibits transitions 
and stabilises single cubes as attractors. As a consequence the coupling acts somehow antiferromagnetic on a coarse 
grained level. 

Performing a coarse graining of the map lattice the resulting symbol or spin dynamics becomes a kinetic Ising model. 
We have been able to identify parameter regions where our dynamical system can be mapped to a finite temperature 



8 



nearest neighbour coupled Ising chain. Depending on the original parameters of the system ferro- or antiferromagnetic 
coupling can be realised, but the ordered phase at zero temperature is always in the antiferromagnetic regime. The 
coarse grained viewpoint also sheds some light on the transient dynamics of the map lattice since the transients 
correspond to a relaxation process in the kinetic Ising model. Therefore, the transient behaviour of the CML is 
related to a non-equilibrium process of statistical physics. 

Within our approach we have successfully linked the dynamics of a coupled map lattice to properties of a kinetic 
Ising model on analytical grounds. Of course our approach is not mathematically rigorous, but we have good indication 
that the results are valid at least in the perturbative regime. The comparison with numerical simulations shows that 
the leading order of perturbation theory is a good description for parameter values e, \S\ < 5 • 10~ 2 . 

For further studies the adaption of the method to coupled maps on a two dimensional lattice seems desirable since 
here also finite temperature phase transitions are possible. This would constitute a further step in the understanding 
of phase transitions in coupled map lattices as exemplified by the Miller Huse model. 



APPENDIX A: 



We want to derive eq. (||) for symbol sequences a., (3 with aj — (3j for j =/= i and ai 7^ fii. First, we remind that the 
single site map f$ is affine on the four intervals (cf . figure 0) : 



K(-2) := [-1, -a], K(-l) := [-a, 0], K(l) := [0, a], K{2) := [a, 1] . 

Consequently, the map T e ^s is affine on the 4^ cuboids 

S 7 := K(rn) x K{ l2 ) ... x K{ 1N ) 
1 ■= {li,l2 ■ ■ -In) with 7j <E {-2, -I, +1, +2} . 

Each cube I a contains 2 N cuboids Sy. The image of a cuboid under T e ,5 is a parallelepiped 

P-y := ^€,6 (5-y) 



(AI) 



(A2) 



(A3) 



which is a weakly deformed cube I a for S-y C I a , because T c= 5 = o (S-y) = I a holds and we are in the perturbative 
regime e, \S\ <§; 1. The distances between the corners of P T and of I a are of the order C(e, S). 
The overlap set OV a} p as defined in eq. (0) then reads 



OV a , = |J (P^nip] 

The intersection of a parallelepiped P-y with Ip can be written as 

N 



(A4) 



P-y n 1/3 = {T e . 5 (x) I x G S-, A H (p 3 [T e , 5 (x)]J = l} 



(A5) 



where 9(x) denotes the Heaviside function. 
Since f3j = aj for j ^ i we have 



= 1 



(A6) 



provided Xj has at least a distance of order 0(e, S) from the endpoints of the interval K (7^). Therefore the set ( |A5| ) 
can be approximated by 



P-yDl p = {T £ , a (x) I x G S y A 9 (A [T e , (5 (x)] i ) = 1 } 



(A7) 



in leading order of perturbation theory. The remaining Heaviside function in eq. ( |A7| ) only depends on the coordinates 
Xi-\, Xi and 2^+1 because of the local coupling of the CML T 6t $. Therefore and since P-y is a weakly deformed cube 
I a , we obtain in the same order of approximation 



P-y n If, = 



{T<3(xW) 

fp(3) 

^ 7.-1 7i 7i+l 



x( 3 ) G 



9(3) a ( B- [T (3) fx( 3 ) > il ^ — 1 \ x r( N - 3 ) 

J 7 I -l7 I 7> + l A yP* I e, 6 V X )\J - 1 J XJ ttitt2...a,-2«i + !-«f 



n J. 



(3) 



x 7(^-3) 



(A8) 
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Here the superscript indicates that the quantities are determined by a map lattice of size N — 3 with x^ 3 ' = 
(ajj-i, X4, Xi-\-i) . The (N — 3) dimensional cube I^a^.-at^ a i+2 ---aN takes the remaining coordinates Xj with j ^ 
{i — 1, i, i + 1} into account. 

If one approximates the map T e ,s by the simplified map 



T e ,«(x) = (1 - e)f s (xi) + ^{f S {xi-i) + fs{x i+ x)) 
T e ,a(x) = /sfo) , V j ^ i , 



(A9) 



one arrives at the result ( A.8) at once. One can use the simplified map T t ,s for the calculation of the pre-image sets 
5 (^K>, p) m leading order perturbation theory. Since "Te,s couples only the coordinates Xi—i, Xi and Xj+i, in this 
approximation also the pre-image sets TJfi(OV a ,p) have the structure of a direct product in eq. (||). Finally, it can 
be easily shown that the blind volume B a of the cube can also be approximated by a direct product of the form 
as in eq. d|). 



APPENDIX B: 



In order to illustrate the main steps for the calculation of the critical values 8 cr it.{e) we focus on the case N ~ 2 

and the transition / ++ — > I |_. Generalisations to N > 2 and different transitions are almost obvious, but require 

some tedious though elementary computations []l2|| . 

For calculating th e overlap set OV_+,++ we introduce the following shorthand notation for the indices of the 
rectangles in eq. (|aJ) 

Si := 5_2i7 S2 := S-22, S3 := S-xi, S4 = S-12 ■ (Bl) 



With the parallelogram Pj := T e< s(Si) the overlap set (A4) reads 



i=l 



Within first order Pi n /++ and P2 H /++ as well as P3 n J++ and P4 n /++ are equal to each other. P3 and the 
intersection P3 n /++ are shown in figure 0. The area of the latter triangular set is e/2 in first order. The intersection 
Pi R i++ is obtained by shifting the just mentioned triangle by an amount —S. 




FIG. 4. The parallelogram P3 = T e ,5(53) and its intersection with the square I++. 

The pre-image set of the overlap, (Oy_+.++), is displayed in figure || where we restrict the parameter range 

to — 2e < 6 < for simplicity. For calculating pre-images of higher order the so called "blind area" B |_ comes into 

play, i. e. the set of points in / |_ which do not have pre-images with respect to the map T ei s. The components of 

(OV r _+.++) in Si and S2 are contained in the subset T of the blind area B [_ (cf. figure |5|). In first order the 

width of the set T is given by the expression 

b T (x 2 ) =f + « 2 , £26 [0,1]. (B3) 
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t Jov )n s 

£,8 - +, ++ 7 



t" (ov )n s, 

e,8 - +, ++ / 1 




T C £ + 



t x (ov )n s 4 

e,8 -+,++' 



T, 8 (oy +i++ )n s 3 



FIG. 5. The pre-image set T7J (OV_ 



blind area -B-+ and its subset T (cf. eq. (B3)) are displayed in gray. 



and its four components in the rectangles Si (— 2e < S < 0). Additionally the 



For convenience in calculating pre-images of higher order we first concentrate on the right rectangles S3 and S4. 
Defining the generations of order k by 



G« := T-}(OV-+,++) D(S 3 US 4 ), 

G (fe) := {x e ( S* 3 U S* 4 ) I T e , 5 (x) e G (fe - 1} } , fe = 2,3, 



(B4) 



figure |(] reveals a beautiful recursive structure of these sets. In order to describe this structure analytically we remark 
that up to first order it is sufficient to compute pre-images with respect to a simplified map (cf. eq. (A9)) 



T £;(5 (x)J = 2x 1 +ef Q (x 2 ) 

T e , 4 (x)l =f (x 2 ) ■ 
J 2 

Then the following properties of the generations G^ , which are inherent in figure [| are easily obtained 



(B5) 



• The first generation G^ consists of two triangles with vertices (0,0), (0,1/2), (— e/2,1/2) respectively (0,1), 
(0,1/2), (-e/2, 1/2). 

• A generation encompasses 2 k triangles each of them having the same area. The area shrinks by a factor 4, 
if one passes from G {k -V to G^ . 

• Two neighbouring triangles of the same generation share a corner or a side with length of order e. 

• The union Eg :— U^ =1 is a simply connected set. 



(k) 

To determine the boundary of X^ we consider its height function 

RW{x 2 ) := inf { x x | (a*, x 2 ) e Y$ } . 
Since i?( fc+1 ) is mapped on R^ by the simplified map T e , we get the representation 

R {k \x2) = -ejZ fJ ^, *2€[0,1] . 
i=l 

For k odd these curves admit 2( k ~ 1 ^ 2 absolute extrema at 

(fc-l)/2 



(B6) 



(B7) 



1 



1 



E 

3=1 



A3 



ij-G {-!,+!}, j = l,2,...(fc-l)/2 



(B8) 



with height 



(fc-l)/2 



E 1 

j=0 



(B9) 
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In the limit k — » oo the set Eg 3 has a fractal boundary, since its construction is analogous to the famous Koch's curve 
J!?]]. The thickness of the set Eg 3 follows easily from eq. (BE) 



MEg) :=sup{M |xGE£} = y 



(BIO) 



1 nG a) /H (i) 
m g <2 7 h (2) 

■ g (3) /h (3) 



-8/2 

FIG. 6. The first three generations and H^ k ' near the right resp. left edge of /_+ . 

A generation £?w has not only pre-images G^ k+1 ^ in the right rectangles S3 and 54, but there are also pre-images 
in the left rectangles Si and S2 (cf. figure ^) 



:= T7j(Oy-+, ++ ) n (SxUS 2 ), 
ff (fc) :={xe (SiUS 2 ) I T £ ,j(x) e G^- 1 )} , = 2, 3, . 



(Bll) 



To re veal the relation between G^ and analytically we just note that for a point y G G^ fc ^ eqs. (B4) and 

flBllD imply that 

T e ,,(x) = T e , 5 (x') = y G G^ , x G G«, x' G . 

Then to first order 

x[ + 1 = -6/2- xx, x' 2 =x 2 (B12) 

follows. Hence, the set is obtained from G^ by a reflection and an additional offset of —6/2 (cf. figure ||). The 
same property follows of course for the limits E^ and Eg 3 . 
If 

Eg 3 C B-+ , (B13) 

holds, no further pre-images of the overlap set OV-+.++ appear, and the union Eg 1 UE^ encompasses all pre-images. 

Consequently, the transition / |_ — > I ++ is not possible, because all pre-images of the overlap set are located near 

the edge of / |_. Therefore, condition (B13) gives the clue for the determination of 5 C rit. ( e )- 

According to eqs. ( B12j ) and ( |B10 ) the thickness of Eg reads 

2e 



ft (Eg) :=su P {1 + .ti |xGEg} = -- 



(B14) 



At the critical value <5 C r;t.(e) one peak at the boundary with maximal height collide s w ith the right border of the set 
T C £_+ (cf. figure @). Since the boundary of the blind area has according to eq. ( |B3| ) a finite slope, the peak with 
the smallest x-i coord inate crosses the right boundary of T at first^J. According to eq. (B8) this peak is located at 
x 2 = 1/3. Then eqs. dml) and (ph yield 



3 This can be shown rigorously with the inequality 



sup{a;i I (xi, xa) £ £j? } < sup {xi \ (xi, 1/3) G Eg 1 } + e ( x 2 - - 
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and consequently we arrive at 



2 £ _ ^(e) = ft (E ^ } = ^ ^ = m = 4 £ (Bi5) 



<W(e) = -^e. (B16) 





J rc 




v 7 




I + 















FIG. 7. Diagrammatic view of the set and the subset T C B-+ of the blind volume for 8 < <5 C i-«.(e). 



In order to show that the transition I |_ — * I ++ is possible for 6 < <5 C rit.(e) the two conditions mentioned at the end 

of section [□] have to be checked. Since \ B |_ is non-empty for 8 < 8 cr it.{e), there exists an open neighbourhood 

of {xi,X2) = (— 1, 1/3) which is contained within the pre-image set T~^° (OV-+.++) for a particular value ko. Next 
pre-images are located near (—1/2, 1/6) and (—1/2, 5/6) and hence enter the inner part of the square I -+. For higher 
generations again four pre-images exist. Therefore it is plausible - and more intricate considerations of |l2| confirm it 
- that the set U^L T~$ {OV-+.++) has a substantial Lebesgue measure. For the second transition criterion we have 
to check whether the points that are mapped into the overlap Oy_+ ; ++ can migrate into the inner part of I ++ under 
further iteration. A point x 6 has a positive x\ coordinate of order 0(e). If one considers the evolution of 

the X\ coordinate (cf. eq. (B5)), its value grows for most points x S OV-+.++ under further iteration, until it reaches 
a value of order one. Therefore, the iterates reach the inner part of /++ after a finite number of steps. 

In conclusion, the transition / |_ — > I ++ becomes possible for 8 < 8 cr u, (e). Computation for other transitions or 

a CML with N = 3 follows the same lines. We stress that the main steps consist in the calculation of images and 
pre-images of overlap sets. The location of the pre-image sets relative to the blind volume determines whether all 
pre-images of the overlap set are located near the edge of the cube only. Consequently, the existence of the blind 
volume influences the numerical value of <5 CT -it.(e) significantly. 



APPENDIX C: 



In this appendix we would like to show that for N — 2 the cubes I |_ and J_| contain attractors in the strict 

mathematical sense, if 8 > — 4e/3. Because of symmetry we can concentrate on the cube I |_. In appendix |B] we 

have shown that the (dominant) transitions / |_ — > I ++ and / |_ — * /__ arc forbidden, as long as 8 > 8 cr u, ~ — 4e/3. 

What remains to be done is that also the off diagonal transition I [_ — > J_| does not appear. 

If 8 < there exists a non-empty overlap set. Considering the four rectangles in / [_ on which the CML T e ,s is 

linear (cf. eq. ([A2])), only the image of the rectangle S-22 '■= [— 1, —a] x [a, 1] intersects the square 7_| for 8 < 0. 
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)-m 




FIG. 8. The overlap set OV-+,+- and its pre-image T~](OV- + ,+-). 
Figure H displays this situation where the parallelogram P-22 = T C! ,5(S'_22) and the overlap set 



OV- 



= P_ 



2 2 



nL 



(Cl) 



are shown. The overlap set has extension ss —6 in the directions of both coordinate axes and hence an area ~ S 2 . 

Note that this area is factor of the order 0(e, S) smaller than those of the overlap sets OV"_+,++ and OV_+ i which 

belong to perturbatively dominant transitions. 

The overlap set (Cl) alone does not ensure for a transition / |_ — > 1^ In fact we show that phase space trajectories 

do not reach this set, so that the transition does not appear. Since we are considering a map with a finite coupling 
e > 0, trajectories do not fill the whole phase space [—1, l] 2 but only the subset T 6t g([— 1, l] 2 ) C [—1, l] 2 . In particular 

points close to the upper left corner of / |_ are not visited. This forbidden domain, previously called the blind volume, 

constitutes the reason why the off diagonal transition does not appear even beyond the perturbation theory. 

To put the argument on a formal level we construct the pre-image sets of the overlap set within the square / |_. 

The first generation set, T~g(OV-+^ ) is located near the corner (—1, 1) of / |_ and has sides of length w —5/2 

(cf. figure ^). As also shown in this figure, the blind area B |_ of the square / y. is also located there and contains 

a square with side length 2e. Therefore, the pre-image set T~j(OV^_+,-| ) is contained in the blind area, as long as 

5 > — 4e. Hence in this parameter regime the overlap set has no pre-image sets T~^(OV_ + _| ) with k > 2, because 

points belonging to the blind area have no pre-images themselves. In particular, the pre-images of the overlap set do 

not intersect the inner part of the square I |_, so that one criterion for the transition / |_ — * J_| is not obeyed and 

the transition is impossible for S > —4 e. 

Summarising, none of the transitions / — (_ — * I a with a 6 { , ++, +— } is possible for S > S cr it. ~ — 4e/3. 

Therefore, in this parameter region attractors in the strict sense reside within the squares I |_ and J_| 
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